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I. INTRODUCTION 



Magneto-transverse light diffusion - more popularly known as the "Photonic Hall 
Effect" (PHE) - was theoretically predicted five years ago by Van Tiggelen Jl], and 
was experimentally confirmed one year later by Rikken This effect is analogous 
to the well-known electronic Hall effect, although the origin of the PHE is somewhat 
different: The PHE finds its origin in the Faraday effect, present inside the dielectric 
scatterers, which slightly changes their scattering amplitude. Most experiments on 
the PHE reported so far used solid samples comprising paramagnetic or diamagnetic 
particles, of size much larger than the wavelength (the Mie scattering regime) 0|. 
The scatterers were embedded in a medium having no magnetic properties, with 
a volume fraction such that multiple light scattering prevailed. A perturbational 
formulation of the Mie scattering of a Faraday-active sphere has been developed to 
explain the origin of the PHE in single scattering |!|. A transport theory for light 
based on this formulation, could produce an estimate of the PHE in multiple light 
scattering, which agreed with experiments Q. Experiments have checked the validity 
of this formulation by measuring the PHE as a function of the volume fraction of the 
scatterers, of the wavelength or of the index of refraction even in the presence of 
absorption 0. 

These experiments can be difficult due to the smallness of the measured PHE. 
It seems therefore natural to try to use other samples such as samples containing 
ferromagnetic particles, in which the magneto-optical effects are expected to be much 
larger than in paramagnetic or diamagnetic samples. The giant magneto-optical 
properties of ferromagnetic compounds were explained by a microscopic theory due 
to Argyres, who derived the form of their polarizability and conductivity tensors jf| . 
This approach was confirmed by Krinchik, whose classical model was able to explain 
both the Faraday effect at optical frequency and the ferromagnetic Hall effect, which 
is the classical Hall effect for the DC magneto-transverse electrical conductivity 

This paper deals with the PHE of liquid or gelled samples of ferrofluids, which are 
colloidal suspensions of ferromagnetic particles. The experiments reported here are 
also the first experimental realization of the PHE in a medium comprising Rayleigh 
scatterers. In these samples, the significant absorption precludes the observation of 
multiple light scattering which does not exclude the PHE. A sketch of the geometry is 
shown in Fig. jjj. This geometry is called magneto-transverse (or Voigt geometry) since 
the magnetic field and the incident light are perpendicular. These two perpendicular 
directions define a plane, and the total scattered intensity above this plane is called 
7 + and the one below this plane it is called I~ . The PHE is a manifestation of an 
anisotropy in the scattered light, which produces a net difference between 7 + and 7~, 
linear in the magnetic field. The quantity r/, which is a measure of the PHE, is defined 
as the difference between these intensities A7(B) = 7 + (B) — 7~(B) normalized with 
respect to the average of the scattered intensities in the absence of any magnetic field 
denoted 7o, 

2(/ + (B)-/-(B)) _ AJ(B) 
' 7+(B = 0)+7"(B = 0) 7 ' 

The difference A7(B) is the photonic equivalent of the difference of potential in 
the electronic Hall effect. In dilute ferrofluids, the ratio r\ can be as high as 10 -3 for 
a magnetic field of lOOOe, which is two or three orders of magnitude above the value 
measured in paramagnetic samples under the same conditions Q. 

This paper is organized as follows: Section [n| recalls some general properties of 



ferrofluids, section III deals with their magneto-optical properties, section [TV] presents 
a theoretical model for the PHE, and section M contains the experimental study of 
the PHE in liquid and gelled samples of ferrofluids. 
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II. FERROFLUIDS 



The stability of magnetic liquids can be destroyed by several factors such as gravity, 
interaction between magnetic moments, gradients in the applied magnetic field, and 
Van der Waals interactions, which favor sedimentation or agglomeration. In practice, 
the stability of the suspension is achieved by chemically grafting charged substituents 
or polymers on the particles, in order to increase the repulsion between the particles, 
of ionic origin (in aqueous ferrofluids) , or of stearic origin (in organic ferrofluids) . 
Because of these stability requirements, the size of the particles in the suspension is 
limited to a radius of a few nanometers H. The ferrofluids used in all experiments 
discussed in this paper, were synthesized by Neveu in the Laboratoire de Physico- 
Chimie Inorganique in Jussieu using a coprecipitation technique first devised by Prof. 
Massart (1981). The samples are suspensions of ferrite particles Fe2CoC>4 dispersed 
in an aqueous ionic solution of citrate. 

In such a suspension, each particle is a magnetic single-domain, of magnetic mo- 
ment 

fi = 4M s fi na 3 /3 » 4, 1 • IOVb, 

where M s ~ 422 kAm" 1 is the saturation magnetization of the bulk material, and 
Hb is Bohr's magneton. From this average magnetic moment and the average dis- 
tance between particles r, a dipolar interaction energy Ed = (/io/47r)/i 2 /r 3 can be 
estimated. For the most concentrated of our samples of volume fraction / = 1%, the 
magnetic interactions are still negligible since at room temperature Ed/ksT ~ 0.17. 
With such a small value, there should be essentially no chain structures in the mag- 
netic liquid. The formation of chain structures in magnetic liquids was predicted 
long ago by De Gennes and Pincus |9L and was studied more recently by Stevens and 
Grest using numerical simulations [lOj. 

In the absence of any magnetic field, the magnetic moments fj, have random orien- 
tation and the net average magnetization of the fluid is zero. When a magnetic field 
is applied, there is a tendency for the particles to align in the direction of the mag- 
netic field. Langevin's theory of paramagnetic gases applies to ferrofluids provided 
that the magnetic interaction are negligible. The behavior of these magnetic fluids is 
known as superparamagnetism g|. The statistical average of the magnetic moments 
~fi is oriented along the magnetic field and its amplitude is 

71 = fiL(u), 

where u = nB /ksT for an applied magnetic field B and L(u) is the Langevin function. 
It is recalled that L(u) — cothu — and L(u) ~ u/3 for u < 1. Quantitative 
comparisons with experiments have shown that a polydispersity function should be 
added to Eq. (^|). Particles sizes obey a log-normal distribution characterized by an 
average radius a and a variance. 

The relaxation of the magnetization is known to have two possible mechanisms: 
either it is caused by a rotation of the particles within the liquid, or by the rotation of 
the magnetic moments inside the particles with respect to their easy magnetization 
axis. The first mechanism is characterized by the Brownian time relaxation H 

Wv 
TB = IcbT' 

where V' is the hydrodynamic volume of the particle, and v is the viscosity of the 
carrier liquid. The second mechanism is characterized by the Neel relaxation time 

1 KV 
t n = — exp( — — ), 

JO Kfli 



where K is the anisotropy constant, V the particle volume, and /o is a characteristic 
frequency of the order of 10 9 Hz. For Cobalt ferrofluids and for a magnetic field of 
frequency ljh, we have tb <s£ uj^ 1 <S tjv. Therefore Neel relaxation is negligible in 
these samples and the dynamic behavior of the fluids follows Brownian relaxation. 
The magnetic moments of the particles in Cobalt ferrite samples remain fixed with 
respect to the particles as confirmed by studies of the linear birefringence in these 
compounds flilj. 



III. MAGNETO-OPTICS OF FERROFLUIDS 

A. Optical properties in single scattering 

Some aspects of the optical properties of ferrite Cobalt samples are now briefly 
discussed. The relative dielectric constant ko in the absence of magnetic field is 
obtained from the dielectric constant of the bulk material, which equals 2.6 — 0.7i 
at Ao = 477nm according to Refs. Jl|[i3]]. The average radius of the particles is 
estimated to be a — 6nm. This corresponds to Rayleigh scattering, since the average 
size parameter, x — 2iran wa ter / — 0.1, which is the ratio of the size of the particle 
to the wavelength in the medium is smaller than one as is also the case for the 
other size parameter y = \fi^xn wa ter ~ 0.2. For these Rayleigh absorbing particles, 
extinction is dominated by absorption because of the behavior of the scattering and 
absorption cross-section as function of x [Q. For x -C 1, and when the imaginary 
part of the index of the scatterer is smaller than its real part (which is the case with 
these ferrofluids), the following relations hold 

/to — 1 8x 4 Kq 



4rfm( — ) > Q sc _ 

Kq + I 6 



KO+2 



The imaginary part of the index can therefore be determined from the absorption 
spectra provided that the real part of the index is also known. Figure ^ repre- 
sents the absorption length as function of the volume fraction on a logarithmic scale. 
The measurements were taken with a spectrometer, and the length of the sample 
was changed so that the absorbance, which is defined as the logarithm of the trans- 
mission, be of order unity. The agreement with the linear behavior predicted from 
Lambert's law, with the slope computed from Rayleigh scattering theory, shows that 
our estimates concerning the index of refraction and the size of the particles (assum- 
ing a monodisperse suspension) are basically correct. The generalization to the case 
of polydisperse samples did not appear to be necessary for this estimate. 



B. Transport properties 

Typical values of the absorption length l a i s of the coherent beam and of the trans- 
port mean free path £* in a sample of volume fraction / = 0.14%, at the wavelength 
Ao = 477nm, are labs = 80/im (measured) and £* = 6cm (theoretically estimated). 
For these samples made of particles smaller than the wavelength, the transport mean 
free path I* is very close to the elastic mean free path I. A treatment of multiple 
light scattering with absorption, becomes necessary when the elastic mean free path 
is smaller than the sample size. For all the samples discussed in this paper, the vol- 
ume fraction is smaller than / = 0.1%, and the condition I 2> L is fulfilled at this 
wavelength, which means that multiple scattering should be negligible. 

Because of the importance of absorption, single scattering should therefore prevail. 
It does not mean however that interference effects are absent, since scattered waves 
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may still interfere as long as they have traveled a distance shorter or comparable to 
the absorption length of the coherent beam £ a bs- This idea will be employed in the 
model of section IVl 



C. Magneto-optical properties 



In the absence of an applied magnetic field, ferrofluids behave as a normal liquid, 
and no birefringence or dichroism, circular or linear, are expected to be present. The 
application of a magnetic field introduces magneto-optical anisotropy, which has been 
the subject of many experimental studies |15lllq]. In the longitudinal configuration, 
in which the wave vector of the incident light is directed along the applied magnetic 
field, the eigenmodes of the electric field are circularly polarized waves. The dif- 
ference in speed and absorption of these waves results in circular birefringence and 
dichroism, which is related to Faraday rotation and ellipticity. The Faraday rotation 
and ellipticity are odd functions of the applied field and are linear at small magnetic 
field as shown in Fig. ^. In this figure, only one curve is plotted, since the curves 
of the Faraday rotation and ellipticity are both proportional to the same Langevin 
function. This was confirmed experimentally by measuring the sample magnetiza- 
tion, the Faraday rotation and ellipticity as a function of the magnetic field at the 
LRCCI laboratory with samples of Fe3C>4. All curves superimpose until volume frac- 
tion of a several percents pju. This experiment confirms that the Faraday rotation 
and ellipticity in ferrofluids have the same physical origin, which is the orientation of 
the magnetic moments in the direction of the magnetic field. 

In the magneto-transverse or Voigt configuration, the eigenmodes of the electric 
field are linearly polarized waves. This implies linear birefringence and dichroism, 
which are both even functions of the magnetic field. Below a critical field of about 
50Oe for the Cobalt ferrite, these effects are quadratic in the applied magnetic field. 
Above this critical field, the effects are linear but still of course even functions of the 
applied magnetic field |[1S| . The linear birefringence and dichroism are also explained 
by Langevin theory pW, as shown by the pioneering experiments of Bacri et al. 
[ pd|pc[ ]. The standard Langevin model however, fails to describe the birefringence 
in spatially ordered samples or in samples in which the particles are not completely 
free to move. This can be the case for instance if the particles are gelled or binded 
to some substrate. The magneto-optical properties of gels were studied theoretically 
and experimentally in magnetic layered silica gels in Refs. Jj^],^. In random sol- 
gels of the type considered in this paper, the four magneto-optical effects have been 
measured 123] . These experiments have shown that the linear magneto-optical effects 
can be inhibited in sol-gels, whereas the circular magneto-optical effects are always 
present and essentially do not change in the gels as compared to the liquids. 



IV. MODEL FOR THE PHE IN FERROFLUIDS 



A. Single scattering T matrix 



A simple model is presented to describe the Photonic Hall Effect in liquid and 
gelled samples of ferrofluids. The samples are dilute and the predominant scattering 
corresponds to single scattering, which is described using a T matrix formulation. 
Our notations are as follows: k and k' are the wave vectors of the incoming and 
outgoing light respectively and B is the direction of the applied field. With each of 
these three vectors, it is possible to associate a local right-handed coordinate system, 
having its z-axis along the wave vector or along the magnetic field. The x and y 
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axis are chosen such that B and k are perpendicular, as is the case in the magneto- 
transverse geometry. 

The incoming beam can be regarded as a linear superposition of plane waves, each 
having a transverse electric field 

Sine = Ei„ c exp [i(kz — tot + 5)] , (5a) 

with 

E inc = Eigi + E 2 g 2 . (5b) 

The wave vector of the light in the medium is denoted k, and the phase 5 can be taken 
as a constant. A convenient choice of polarization vectors is: gi = g = k x k'/|k x k' | , 
g2 = k x gi. The scattered beam can be decomposed in the polarization vectors 
g 2 = k' x gi and gi = gi. This choice of polarization vectors was introduced by 
Van De Hulst: one polarization vector is located in the scattering plane while the 
other one is perpendicular to it |Q. This choice is not well defined in the particular 
cases of forward and backward scattering because there is no scattering plane in these 
cases, but these two particular configurations are not relevant for the PHE. In the 
frame of the particle, in which quantities are denoted with a prime, the scattering 
matrix reads 

/ t it'x \ 

(6) 

The complex- valued coefficients t' , t\ and t' 2 describe a point-like scatterer in a 
magnetic field. 

To show the relationship between the properties of a single particle and the macro- 
scopic magneto-optical properties of an ensemble of particles, a statistical average of 
the magnetic moments of the particles has to be performed, as discussed in Ref. 
[ p6[ . When the coupling between the magnetic moments is neglected, and when a 
random distribution of scatterers is assumed, the average T matrix has the same 
form as Eq. ([]) with the following parameters to — t' + t' 2 L(u)/u, t\ = t[L(u) and 
t 2 — t 2 L 2 (u) = t' 2 (l — 3L(u)/u), where u — fiB/ksT for an applied magnetic field B, 
since 

/ t' +t' 2 L(u)/u it[L(u) \ 

t=F={ -itiL(u) t' + t' 2 L(u)/u . (7) 

V t' + t' 2 (l-2L(u)/u) J 

This result implies that the Faraday rotation associated with ti, is proportional to 
the Langevin function L(u), and that the linear birefringence is associated with the 
second Langevin function L 2 (u) contained in t 2 , as confirmed by experiments. 

The scattering matrix, which relates linearly the outgoing and the incoming electric 
field needs now to be expressed in the basis of the polarization vectors introduced 
above. This yields to the following result: 

Su = t + t 2 (B ■ g) 2 , (8a) 
(B • k') (B ■ g) 

S12 = -t 2 ± ; M Sj , (8b) 

Ik X k'l 



Ik x k' 
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. „, (k-k')(B-k') . . 

S22 = tok-k' + t 2 - — — ^--itiB-tkxk'). (8d) 

|k x k' 

These expressions are valid only for the magneto-transverse geometry but at any 
order in the magnetic field. They can be simplified by introducing a coordinate 
system in which the wave vector of the scattered wave k' makes an angle 8 in the 
scattering plane with respect to the incoming direction and an azimuthal angle ip 
with respect to B as shown in figure Q With this choice, 

/ sin 6 cos ip \ 
k' = I sin sin ip I . 
V cos 9 J 



B. PHE for unpolarized incident light 

The scattering matrix of Eqs. (p^-pd|) describes single scattering. This approach 
is not sufficient to account for the PHE in ferrofluids, because the PHE of a single 
magnetic Rayleigh scatterer is zero as was found in Ref. |3|. It is therefore necessary 
to take into account interferences among different single scattering events, in order 
to explain the observed non-zero value of the PHE. To this end, one can start with 
the scattering matrix of a set of particles at position r; and of moment fn, which can 
be written 

(S) = 5^S( W )|r i )<r i |. (9) 

i 

The summation can be performed with the following approximation. Because of 
the significant absorption in the samples, it is meaningful to neglect all interferences 
effects outside a range of the order of the absorption length of the coherent beam 
tabs- In a medium of size l a bs, which is assumed to be spherical, interferences can be 
taken into account exactly in single scattering since I* 3> l a b a - In the basis of Van 
De Hulst for the polarization, the scattering matrix has the following form 

{S ij )»nJ ab \^- k '^S ij d 3 r i = f^f-yG(v)S ih (10) 

where n is the number of particles per unit volume, G{v) = 3(sin v — v cos v) /v 3 is the 
phase function for Rayleigh-Gans scattering and v = 2k£ a t, s sin(#/2). Relation ( px| ) 
expresses a well known result in optics: the proportionality of the scattering matrix 
of a Rayleigh-Gans to the scattering matrix of a Rayleigh scatterer |24|. From the 
scattering matrix, the scattering cross-section can be computed for the given states 
of incident and outgoing polarizations. The scattering cross-section, averaged over 
incident and outgoing polarization, reads: 



»,J=1>2 



To compute the PHE 77 defined in Eq. (|l|), it is necessary to integrate the scattering 
cross-section with respect to the direction k' including a magneto-transverse projec- 
tion factor B ■ (k x k'). This factor allows to calculate the light flux projected onto 
the magneto-transverse direction B x k, which is the direction of the normal of the 
detector in the experiment. After this integration, only terms odd in the magnetic 
field remain. For unpolarized incident light, the PHE finally reads 



6 



with 



Q'm(ioti) 
71 = ^^' 

The integrals in Eq. (|l^) can be calculated in the limit k£ a bs 1, which gives 

_ 371 \n(ki abs ) 

Vunpol — ~ , „ 

It is interesting to note that when the size parameter of the particles goes to zero 
x — *• 0, the parameters 71 and r] unpo i are independent of the size of the particles. This 
makes this limit insensitive to polydispersity, which is always difficult to estimate in 
this kind of experiment. 



C. PHE for polarized incident light 



Polarized incident light may be either circular or linear. A state of linear polariza- 
tion is described by 

Ei = cos/3, E2 = sin/3, 

in Eq. (|5b[). A sketch of the geometry is shown for this case in Figure ^. When (3 = 0, 
the electric field of the incident light is parallel to the applied magnetic field. 

Contrary to the case of unpolarized PHE, several terms present in the scattering 
matrix now contribute. After having carried out the integration over the outgoing 
wave vector and the average with respect to the outgoing polarization, the numerator 
A/(B) and the denominator Iq of the PHE take the form 



AJ(B) = tyA 



--Im(t 2 tt) cos(2/3) + 2Im(t tt) sin f3 



Io = 2A\t \ 



d6 sin 3 6 cos 9G 2 [v{6)\ , 
d6 [cos 2 13 + cos 2 9 sin 2 0\ sin 2 9 G 2 [v{9)\ , 



where A is a known constant of proportionality. In these expressions, the integrals 
over the scattering angle 9 can be done in the limit k£ a t s 2> 1. For the last equation, 
this gives 



Io = 2A\t \ 



d9 sin 2 9 G 2 [v(0)] , 



which means that the denominator Io of the PHE is independent of the angle (3 char- 
acterizing the state of linear incident polarization. The PHE for polarized incident 
light now reads: 



r](f3) = 2r] unpo i 



sin 2 (3- ^- cos(2/3) 
471 



where r) unpo i is the value of the PHE obtained from Eq. (114) f or unpolarized incident 
light, and 72 has been defined in the same way as in Eq. (p-3[) by 



72 



Im(t 2 t? 

M 2 



Averaging over (5 in Eq. ( JlTj ) reproduces the value fjunpol- Two particular cases of 
Eq. ( |17[ ) will be of interest: when the linear magnetic birefringence is negligible, one 
has 

r){P) = 2r/ unpo i sin 2 (3. 

On the other hand, for a dominant contribution of the birefringence, one gets 

,a\ ^ 2 oa\ 372 ln(fc4i, s ) 
V{P) = -Vun P oi—cos{2f3) = \- '- cos(2/3). 

The dependence of the PHE on the state of circular polarization can be done very 
similarly. A state of circular polarization can be written 

F - 1 F - ±l 
till — —7=, £ji — — 7=- 

One finds that the PHE is independent of the state of circular polarization left or 
right of the incident light 

This value corresponds to the case of linear polarization at the angle j3 = 7r/4, or to 
the case of unpolarized incident light. 



V. EXPERIMENTAL RESULTS 



A. Dependence of the PHE as function of the field and concentration 

The light intensities I + and I~ , represented in Fig. ^ travel trough optical fibers 
from the sample to the detectors which are photo-diodes or photo-multipliers. The 
difference A7(B) is measured using a lock-in. More details on the experimental setup 
can be found in Ref. Q. The measurement of the phase difference of AJ(B) with 
respect to the magnetic field defines the sign of the PHE. This phase for a true signal 
should be therefore either or ir radians. The sign of the PHE is a very important 
feature which has been seen to be very sensitive to several experimental parameters 
(Verdet constant, index of refraction, concentration, polarization..). 

Our first experimental study addresses the dependence of the PHE as a function 
of the applied magnetic field. As shown in Fig. Bl the PHE is linear in the applied 
magnetic field at low field. In ferrofluids, the deviation away from a linear behavior 
sets in at fields higher than a few hundreds of Oe. Experiments done at higher field 
show the same saturation which is observed in figure H for the Faraday effect. This 
aggrees with the prediction of Eq. dlif ) and Eq. (^), and confirms that the PHE 
is produced by the orientation of the magnetic moments according to the Langevin 
model. Unless specified otherwise, the rest of the measurements of the PHE was done 
in the linear regime for the PHE. 

The variation of the PHE with volume fraction /, is shown in Fig. ^ for a gelled 
sample, which was prepared with about 1% of gelatine in volume. The experiment 
was repeated with liquid samples of ferrofluids. In both cases, a linear behavior with 
respect to the volume fraction was found. The dispersion of the experimental points 
in Fig. ^| around the linear behavior is believed to be due to the evaporation of water 
in the gelled samples, and was absent in the experiment with liquid samples. The 
approximation used in Eq. (jl4j) is valid since k£ a t, s — 1.6T0 3 for the most concentrated 
of the samples used of volume fraction / = 0.14%. Since the absorption length £ a i, s 
is inversely proportional to the volume fraction of the scatterers, Eq. (FL4J) predicts a 



8 



linear dependence of r\ with respect to the volume fraction /. This must be caused 
by a difference AJ(B) quadratic in f 2 while the scattered intensity Jo is proportional 
to /. A quadratic behavior for f 2 in AJ(B) means that the PHE involves more than 
one particle, because the probability of scattering for a photon is proportional to /. 
These experiments confirm our previous statement that a single Rayleigh scatterer 
is not able to generate the PHE but that more than one Rayleigh scatterers or one 
scatterer of finite size is necessary to get a non-zero value of the PHE [0 . 

The PHE per unit magnetic field in the experiment of Fig. ^ is estimated to be 
(Fri/dBdf = 3.1tT 3 T -1 . In order to compare this value with a theoretical estimate, 
it is necessary to evaluate the coefficients of the scattering matrix. For a low applied 
magnetic field, these coefficients can be determined from the anisotropy in the index 
of refraction [E5j. The dielectric constant of the sample is assumed to have the 
following form at low applied magnetic field: 

£ij = toSij + ieFtijkBk- (22) 

Using our measurements of the Faraday rotation and ellipticity in this sample, we 
estimate the antisymmetric part of the dielectric constant ef to be cLef/cLE ~ (—0.5+ 
1.5i) • 10" 4 T~ 1 . With this value, one can calculate the antisymmetric part of the 
relative dielectric constant kf of a single particle. The parameter 71 can then be 
obtained from the relation 



Im(totf) 

71 = n*r~ = ~ 3Im 



Hp 



(k - 1)(«o + 2) 



This yields d 2 rj/dBdf ~ 1,2 • 10 -3 T -1 , which is to be compared with the slope 
of Fig. ^ d 2 r)/dBdf ~ 3 • 10 -3 T -1 . Because of experimental uncertainties due to 
the lack of stability of the gel, this can be considered as a fair agreement. The other 
experiments with gels which are reported in this paper, were done with sol-gels which 
are stable. 



B. Experimental study of the polarization of the PHE 

The dependence of the PHE on the volume fraction is an important test of the 
model. However a lot of physical and optical properties depend on the volume frac- 
tion, which makes the interpretation of this kind of experiments difficult. The polar- 
ization is a much more robust quantity to consider, since it should be independent of 
many experimental conditions. This section adresses the dependence of the PHE on 
the incident polarization for liquid and gelled samples of ferrofluids. The polarization 
of the scattered light is still averaged by the detector, which is not sensitive to the 
polarization of the light, as was checked experimentally. The opposite experiment, 
consisting of sending unpolarized light and measuring the polarization of the out- 
coming light, is equally interesting, but has not been carried out, because it is more 
difficult to realize. 

The first experiment deals with sol-gels of ferrofluids. The gels were produced in the 
group of J. A. Serughetti from Departement de physique des Materiaux de l'Universite 
Claude Bernard in Villeurbanne. The gels were made from the isotropic liquid state 
with no magnetic field present during the gelation. During the measurement of the 
PHE, a magnetic field of 250Oe and of frequency 560Hz was applied. In Fig. ^a 
the normalized PHE is shown as a function of f3, together with the prediction from 
Eq. (|li|), when magnetic linear birefringence, the second term in Eq. (1^), is assumed 
to be absent. This assumption is consistant with the observed weak sensitivity of 
this curve with respect to the amplitude of the applied magnetic field. We have 
checked that the dependence of the PHE on /3 for linearly polarized incident light 
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comes in only from the numerator of r\ in Eq. ([fj), as no dependence was found in 
the incoherent background Iq as predicted by the model of the previous section. 
The magnetic birefringence was assumed to be absent in the comparison with our 
experiments since the it should be inhibited in the gel. Although the general features 
of the curve are captured by the model, there remains a clear discrepancy between 
experiments and theory. This discrepancy follows a cos f3 2 sin 2 (3 law, which probably 
indicates that this point might be explained by taking into account higher scattering 
orders, such as double scattering for instance. 

Figure 0b shows measurements of the PHE as function of the angle j3 together with 
the theoretical prediction of Eq. (§o|). The sample is a liquid Cobalt ferrofluid, of 
volume fraction 0.025%, and the applied magnetic field is 40Oe. It can be first noted 
that the observed dependence of the PHE is completely different from the previous 
case: the symmetry of the curve with respect to the zero axis is different. There 
is now a range of values of (3 for which the PHE changes sign, and the effect no 
longer vanishes at /3 = 0, but is maximum. Contrary to the case of gelled samples, 
the experimental result is now compared to the theoretical model for dominant linear 
magnetic birefringence. In these ferrofiuids, the particles are indeed free to reorganize 
themselves and acquire shape anisotropy, which results in a linear magnetic birefrin- 
gence generally larger than the Faraday effect. In the same experimental conditions 
of figure (^b) and in a sample length of 3mm, the Faraday rotation was measured to 
be 0.1°, to be compared to 1.7° for the linear magnetic birefringence. The magnetic 
linear birefringence corresponds to the term 72 in Eq. (po[), which is an odd function 
of the magnetic field very much like 71 in the unpolarized case. At low magnetic 
field, 72 is proportional to B 3 , and should therefore generate an oscillation at the 
modulation frequency of the magnetic field. This time, a very good agreement with 
the predicted law in — cos (2/3) is found. 

The experiment was repeated with circularly polarized incident light, in liquids and 
in gels, and fully confirmed the result of Eq. (pl|). There is virtually no measurable 
difference between 77 + et r\~ . The value of r/ + corresponds indeed to the value of the 
PHE with linear polarization for an angle /3 = 7r/4. 



VI. DYNAMICAL ASPECTS OF THE PHE 



In the last section of our study of the PHE in ferrofiuids, some dynamical aspects 
are discussed, in relation to the rotation of the magnetic moments of the particles 
induced by the magnetic field in ferrofiuids. In section [n], it was noted that the 
main process of relaxation of the magnetic moments in the Cobalt samples is the 
brownian mechanism. Since the Faraday effect is produced by the orientation of the 
magnetic moments, it should follow the same relaxation law as the magnetization. 
This relaxation of the Faraday rotation as a function of the frequency of the magnetic 
field was indeed observed |27]]. The characteristic relaxation time can be obtained 
either from the study of the frequency dependence of the Faraday rotation or from 
the study of the transient response as the magnetic field is turned off. Both methods 
provide comparable values of the relaxation time. From Eq. (Q), it can be noted that 
the relaxation time for brownian relaxation is directly linked with the viscosity of 
the suspension. Using this idea, Bacri et al. in Paris and B. Payet in St-Etienne 
have experimentally shown that measurements of relaxation times could be used 
to determine the viscosity of the carrier liquid v. A viscosimeter based on this 
principle could achieve an accuracy of a few percents and a range of several decades of 
viscosity. More recently, measurements of relaxation times are now used for biological 
applications of ferrofiuids [ p8[ . 

The magnetic susceptibility of ferrofiuids is well described by Debye theory, orig- 
inally applied to dielectrics but also very useful for magnetic liquids. Debye theory 
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essentially relies on the following hypotheses: 



All magnetic moments of the particles are fixed with respect to the particles 
All particles have identical dimensions and magnetic moments 
There is no interaction between the particles. 

Within Debye's model, the magnetic susceptibility is 

X("h) = 1 *° = \x(u H T)\e*^\ 
1 + iujht 

where luh stands for the frequency of the magnetic field and r for the relaxation time 
of the liquid. The proportionality of the magnetic susceptibility ( p3| ) and the Faraday 
rotation (measured both in amplitude and phase) has been established in Ref. j27[ . 
Therefore, it is expected that the PHE should be also proportional to this complex 
magnetic susceptibility. It is very important not to confuse the phase ^(luht) defined 
in Eq. ( |23| ) with the phase of the scattered light, which is often written as a complex- 
valued quantity. The phase </?(o)£rr) is defined unambiguously with respect to the 
magnetic field and varies at the frequency of the magnetic field uih, whereas the 
phase of the scattered light oscillates at optical frequency. 

In order to show the applicability of Debye's theory to the PHE in ferrofluids, the 
amplitude and phase of the PHE were measured as function of the product luht, as 
shown in Fig. ^. Experiments were carried out in samples of different viscosities by 
diluting a suspension of (liquid) ferrofluid with water, which was originally a solution 
in glycerol. Three samples of identical volume fraction / = 0.025% in ferrofluid were 
used: the original solution with 100% of glycerol of viscosity 1093Pa.s, a solution with 
90% of glycerol of viscosity 384Pa.s, and a solution with 75% of glycerol of viscosity 
60Pa.s at 25°. A magnetic field of lOOOe was applied during the experiment. The 
curves for different viscosities do superimpose on each other, which confirms that the 
PHE only depends on the product loht. 

The saturation of the PHE as function of the magnetic field and the measurement 
of relaxation time both confirm that the origin of the PHE resides in the orientation of 
the magnetic moments of the particles, just like for the Faraday effect. The observed 
behavior as function of frequency in figure [s] is compared with the prediction of the 
Debye theory of Eq. (p3j). The measurements of the amplitude and phase of the 
PHE agrees with the Debye model when the parameter lobv is below a value of 
about 10 5 Pa. Above this value, a deviation is seen between the experiment and the 
theoretical curve of the amplitude of the PHE, and a more dramatic deviation is seen 
in the phase at high frequency. The experimental setup should not be responsible 
for this deviation, even at the highest frequency of lob = 1.875kHz. The PHE at this 
frequency is about two orders of magnitude smaller as compared to its value at low 
frequency, but the noise level is still very low for these measurements (less than 1%) . 
In these rather dilute samples, the origin of this deviation, the caracteristic time of 
which is about 8ms, is unclear. 



VII. CONCLUSION 



We have presented a theoretical model of the PHE, which captures the essential 
features of the experiments reported in this paper. This work confirms our under- 
standing of the PHE in a system of Rayleigh absorbing scatterers, and its connection 
with the Faraday effect. In ferrofluids, as shown by static and dynamic measurements, 
the Faraday effect and the PHE were both found to be produced by the orientation 
of the magnetic moments. This work on the PHE in ferrofluids could also be useful 
to improve our understanding of complex media containing ferrofluid particles, such 
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as media in which the particles are not free or are only partially free to move. The 
polarization dependence of the PHE was found to be sensitive to the degrees of free- 
dom of the particles in the medium. In the magneto-transverse geometry, this type 
of experiment could provide information complementary to the one provided by the 
magnetic linear birefringence. 

We wish to thank H. Roux, C. Bovier and G. Rikkcn for their help in the experi- 
mental part, J. Ferre, R. Perzynsky, and J-C. Bacri for the enlightening discussions 
on ferrofluids, and T. Lubensky for stimulating discussions. 
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FIG. 1. Sketch of the geometry of the PHE. The sample is shown as a cylinder. The scattered intensi- 
ties above (and below) the plane containing the incident light direction and the magnetic field are called 
respectively I + and I~ . The difference AJ(B) is the photonic equivalent of a difference of potential in the 
electronic Hall effect. 
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FIG. 2. Logarithm of the absorption length logio(£ a 6s/l/im) as function of the logarithm of the volume 
fraction / at the wavelength Ao = 477nm, for several aqueous solutions of Fe2Co04. Points represent 
measurements and the line results from Lambert-Beer's law (the slope is estimated from Rayleigh scattering 
theory for a monodisperse distribution of particles of radius a = 6nm using the index of refraction of the 
material at this wavelength 2.6 — 0.7i.) 
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FIG. 3. Normalized measurements of Faraday rotation and ellipticity in a sol-gel of Cobalt ferrite of 
volume fraction 0.04%. The superposition of the normalized curves for the two effects is an indication 
of their common origin, both effects scaling with the Langevin function introduced in section For the 
maximum field 82240e, the Faraday rotation is 0.7° and the ellipticity 1.91° in a sample holder of length 
2mm. The wavelength in this experiment was Ao = 633nm, and these measurements were taken in a static 
magnetic field. 





FIG. 4. Sketch of the geometry in the frame of the laboratory. The vector B defines the axis x, the 
magneto-transverse direction is along the axis y, and the incident wave vector is along z. The angle j3 
characterizes the state of linear polarization of the incident electric field. 
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FIG. 5. (a) The squares represent measurements of the PHE 77 at low field, as function of the applied 
magnetic field. The curve only connects the experimental points. It can be noted that the PHE is linear in 
magnetic fields smaller than about a hundred of Oe. This experiment was done with the same sample than 
in Fig. (0). (b) Same experiment with similar samples but with higher magnetic field, so that the saturation 
of the PHE is now clearly visible. The circles are experimental points and the line is a fit using the Langevin 
function. 
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FIG. 6. Points are experimental measurements of PHE r\ as function of the volume fraction /. Samples 
are gels made with 1% of gelatine in volume and Cobalt ferrofluid. A linear fit gives the slope per unit of 
magnetic field cPrj/dBdf = 3 • 10 _3 T _1 . 
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FIG. 7. (a) Normalized PHE for a sol-gel sample of ferrofluid as function of the a ng le f3 of the incident 
linear polarization. Points are experimental and the curve is a prediction from Eq. (M). A magnetic field 
of 250Oe and frequency 560Hz was applied on the same sample studied in Fig. Ha. (b) Normalized PHE for 
liquid ferrofluid as function of /3, together with the prediction of Eq. (^). A magnetic field of 40Oe and 
frequency 20Hz, was applied on a sample of volume fraction 0.025% of aqueous Cobalt ferrite ferrofluid. 
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FIG. 8. Amplitude (a) and phase (b) of normalized PHE as a function of the product of the frequency of 
the magnetic field u>b and the viscosity of the suspension v for samples of different viscosity: the original 
solution with 100% de glycerol (squares), a solution with 90% of glycerol (circles), and a solution with 75% 
of glycerol (diamonds). The solid curves in (a) and (b) are fits using Eq. (j%|). The experiment has been 
realized with a field of lOOOe and with samples of volume fraction 0.025%. The Debye model predicts a 
behavior of the phase at high frequency, which is not seen in the experiment. 
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